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ABSTRACT 

With this paper we introduce the concept of apparent structure of a GRB jet, as opposed to its 
intrinsic structure. The latter is customarily defined specifying the functions e{9) (the energy 
emitted per jet unit solid angle) and r(0) (the Lorentz factor of the emitting material); the 
apparent structure is instead defined by us as the isotropic equivalent energy Eiso{dv) as a 
function of the viewing angle 9^. We show how to predict the apparent structure of a jet given 
its intrinsic structure. We hnd that a Gaussian intrinsic structure yields a power law apparent 
structure: this opens a new viewpoint on the Gaussian (which can be understood as a proxy for 
a realistic narrow, well collimated jet structure) as a possible candidate for a quasi-universal 
GRB jet structure. We show that such a model (a) is consistent with recent constraints on the 
observed luminosity function of GRBs; (b) implies fewer orphan afterglows with respect to 
the standard uniform model; (c) can break out the progenitor star (in the collapsar scenario) 
without wasting an unreasonable amount of energy; (d) is compatible with the explanation of 
the Amati correlation as a viewing angle effect; (e) can be very standard in energy content, 
and still yield a very wide range of observed isotropic equivalent energies. 

Key words: radiation mechanisms: non-thermal - relativistic processes - gamma-ray burst: 
general 


1 INTRODUCTION 


In a 1999 preprint Lipunov, Postnov and Prokhorov introduced, 
possibly for the first time, the idea that Gamma-Ray Bursts (GRBs) 
could be “standard energy explosions”, i. e. events with a standard 
energy reservoir. In t heir final article, published two years later 
( [Lipiinov et al.|200T l, the authors identified i?o ~ 5 x 10®^ erg 
as a plausible value for this universal energy, and they described 
two possible scenarios: in the first, the standard energy is emitted 
inside a conical jet whose semiaperture 6'jet varies from one GRB 
to another (Fig.[T]a); in the second, the beam pattern, made up of 
two coaxial conical components and an isotropic component, is the 
same for all GRBs (Fig. [T]b). In their view, the wide range of ob¬ 
served luminosities of GRBs could be accommodated in either the 
first scenario, with the brightest events being the most collimated, 
or the second picture, with the viewing angle being crucial to de¬ 
termine which part of the beam mostly contributes to the observed 
fluence. 

Soon later |Frail et al.| ( [2001| > analyzed a sample of 17 Long 
GRBs for which a jet break in the afterglow light curve was 



Figure 1. Sketch of the two possible scenaiios described in |Lipunov et al.| 
j2001j . In the first one (a) GRB jets are seen on-axis, and they differ by 
their semiaperture and consequently their observed energy, with the bright¬ 
est being the most collimated; in the second one (b) the jet configuration is 
such that the viewing angle determines which component contributes most 
significantly to the received energy. 
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identified, and thus a measure of 6'jet was available, finding a 
surprising clustering of the collimation-corrected energy E-f = 
Eiso (1 COS 0jet) « i5iao0fet/2 around the universal value E-/ ~ 
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5 X 10““ erg (which implies a correlation Eiso oc ). The re¬ 
sult was interpreted as evidence that the emission is in fact beamed 
inside a conical jet: this supported the first scenario proposed by 
|Lipunov et al.| ( [2001| l, finally tracing the very wide range of ob¬ 
served isotropic equivalent energies of GRBs to a single “real” 
value. 

Next year, [Rossi et al.| ( [2002| l and [Zhang & MeszarosH2002) 
interpreted the same result in a different way, closer to the second 
scenario proposed by [Lipunov et al.| ( [2001) : their claim was that 
the correlation Eiso oc was instead due to the existence of a 
universal jet structure, with the jet energy per unit solid angle being 
£(6^) oc 9~^. This particular energy configuration, along with the 
assumption of a strong relativistic beaming of the emitted radiation, 
implies that Eiso oc 6~^, where 9^ is the viewing angle. Based 
on a simulation o f the afterglow ligh t curve^roduced by such a 
structured jet (SJ) Rossi et al. 1 2002 1 claimetjjthat 9^ ~ Sjet, and 
thus Eiso oc 9^^. 

The idea of a quasi-universal jet structure (that is, univer¬ 
sal with some dispersion of the structure parameters) stimulated 
a number of papers exploring the idea and its consequences. Here 
are some examples: [Granot & Kumar|j2003[ > and [Kumar & Granot[ 
( [2003^ constrained the possible jet structures by a qualitative com¬ 
parison of simulated afterglow light curves with the observed ones; 
[Rossi et^T] ( [2004 ^ showed that polarization measures could be a 
crucial tool to discern between the SJ and the uniform jet; [Zhang[ 
[etaLjpOO?) and |Dai & Zhang[ ( |2005[ l proposed the interpretation of 
X-Ray Flashes (XRFs) and GRBs as a unique population of bursts 
with a Gaussian structured jet, and tested this hypothesis against 
many observational constraints, finding general consistency; on the 
contrary, |Lmnb^e£^ 1 2005 [> fou nd that the universal structured jet 
proposed by Rossi et al. ( 2002\ fails to predict the right observed 


number ratio of XRFs to GRBs; relativistic hydrodynamical sim¬ 
ulations (e.g. [Zhang et al.[j2004[ >; [Morsony et al.[j2007^ ) showed 
that the interaction of the jet with the stellar envelope prior to the 
break out (in the collapsar scenario of Long GRBs) leads inevitably 
to some structure in the jet properties, but it remains unclear if this 
structure is likely to have any degree of universality;[Lundman et al.[ 
( [2013^ , within a photospheric emission model, obtained a low en¬ 
ergy photon index a consistent with the observations assuming a 
SJ where the Lorentz factor r(0) varies as a power law with the 
angular distance from the jet axis; we ( [Pescalli et al.[2Q15t recently 
showed that the observed luminosity function of Long GRBs is 
consistent with the SJ model, provided that the energy structure 
is much steeper than the original 9~'^ (at least 9~'^ seems to be 
necessary). 

Despite all these efforts and some successes, no consensus 
about the viability of the quasi-universal SJ hypothesis has been 
achieved so far. 


predictions of the SJ model with the observations. In the follow¬ 
ing sections, we will try to make the distinction between intrinsic 
and apparent structure as clear and as rigorous as possible; for the 
moment, suffice it to say that the intrinsic structure here is under¬ 
stood as the energy emitted by different portions of the jet at differ¬ 
ent angular distances 9 from the jet axis, while the apparent struc¬ 
ture describes the energy received by observers that see the same 
jet under different viewing angles 9^. While the intrinsic structure 
can be due to the jet formation process itself (e. g. the [Blandford] 
[& Znaje^ ( [1977^ process) or to the subsequent interaction of the 
jet with the stellar envelope (in the collapsar scenario), the appar¬ 
ent structure depends on how relativistic beaming effects shape the 
emission from each part of the jet. From an observational point of 
view, it is the apparent structure that determines the isotropic equiv¬ 
alent energy, the observed luminosity function and the like; from a 
theoretical point of view, one would like to reconstruct the intrin¬ 
sic structure to find out e. g. the actual energy emitted by the jet 
and, through the efficiency, the total energy (kinetic plus internal 
and possibly magnetic) of the jet itself; thus, a clear distinction be¬ 
tween the two, and some insight on their interdependence, are to be 
worked out. 

The importance of such a distinction was partly pointed out 
in a remarkable work by [Graziani et al.[j2006[ l, but their study as¬ 
sumed a constant bulk Lorentz factor profile r(0) = F. Few works 
to date (as far as we know) assign a variable Lorentz factor pro¬ 
file r(0) to the jet (e. g. Kumar & Granot 1 2003^; Lundman et al. 


( [2013^ ) and none reports predictions about the apparent structure 
of a GRB jet within such a model. 


The structure of the article is the following: 

(i) in !|^and |3.1[ we give a rigorous definition of intrinsic and 
apparent structure, and we make some examples to clarify the con¬ 
cepts; 

(ii) in §3.2[ we introduce two formulas to compute the apparent 
structure and the spectrum of a SJ given its intrinsic structure; in the 
following subsection, we compare the predictions of our formulas 
with previous treatments and show that they are consistent; 

(iii) in what follows next, we analyze the particular case of a 
Gaussian intrinsic structure, showing that (under the assumption 
that also the Lorentz factor has a Gaussian profile) its apparent 
structure is not Gaussian, but rather it is quite well described by a 
power law; we then show that a Gaussian quasi-universal jet model, 
with very reasonable parameter values, is consistent with recent 
constraints from the observed luminosity function of GRBs; 

(iv) we show that the model is consistent with the Amati corre¬ 
lation being a viewing angle effect; 

(v) in the Appendix we give detailed derivations of the formulas 
presented in §3.2[ 


1.1 Aim and structure of this paper 

The aim of this article is to introduce the concept of apparent struc¬ 
ture and to show that it is a necessary tool to correctly compare the 


^ In [Rossi et al. ]S, the afterglow light curve of their SJ seen under a 
viewing angle 6-v is found to show a feature similar to the jet break predicted 
for the uniform jet, with the coincidence that the break time is approxi¬ 
mately the same as that of a uniform jet, seen on-axis, with semiaperture 
^jet = The correspondence is not exact, as discussed two years later in 
[Rossi et S^j2004|, but the difference should be small in most cases. 


2 INTRINSIC STRUCTURE 

Following [Rossi et al.[j2002} and [Zhang & Meszaros[ p002^ we 
define the intrinsic structure of the jet as follows: 

• we set up a spherical coordinate system with the central engine 
at its origin and the jet directed along the z axis; 

• we define the function t{9) as the energy emitted (during the 
prompt emission) by the portion of the jet comprised between 9 
and 9 + d9, divided by the corresponding solid angle, i. e. t{9) = 
Tj dE(9)/2Tt sin 9 d9, where dE here stands for the total energy 
(kinetic plus internal and possibly that of the magnetic field) of the 
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Figure 2. Example apparent structure of a uniform jet in the ultrarelativistic 
limit (black dashed line) and for F = 200 (red solid line). The isotropic 
equivalent energy is normalized to its maximum value, con'esponding to 
the jet observed on-axis. A jet semiaperture 6jet = 3° is assumed. 


jet portion, and rj is the prompt emission efficiency, which might as 
well depend on 

• we assign a Lorentz factor r(0) to the emitting material com¬ 
prised between 0 and 6 + d9 during the prompt emission. 

The functions e{9) and T{6) then define what we call the intrinsic 
structure of the jet. 


3 APPARENT STRUCTURE 


3.1 Definition 


We introduce here our notion of apparent structure. Let 6^ be the 
viewing angle of an observer looking at a GRB jet, i. e. the angle 
between the jet axis and the line of sight. We call “apparent struc¬ 
ture” the function Eiao{dv), namely the isotropic equivalent energy 
inferred by the observer as a function of 8^. For the sake of clarity, 
let us apply this definition to some examples: 

(i) an isotropic explosion, defined by £(0) = e G [OjTr], 
would clearly have 

Eiso = 47re (1) 


for all viewing angles; 

(ii) the “classical” uniform (“top-hat”) GRB jet has 


6 ( 0 ) = 


and 

rw = I 


e 8 < 6jet 

0 8^ 8jat 


T 8< 

1 8 ^ 0jet 


( 2 ) 

(3) 


In the ultrarelativistic limit (E —>■ cxd) the uniform jet is indistin¬ 
guishable from an isotropic explosion as long as 8v < 8jet ,because 
the relativistic beaming prevents ( |Rhoads|1997| l the observer from 
“seeing” anything not expanding exactly along the line of sighQ 


^ The implicit assumption is that the jet expansion is purely radial with 
respect to the central engine. 


For the same reason, the GRB is always undetected if 8^ > 9jet - In 
other words, the apparent structure (dashed black line in Fig.|^ is 


EiBo{8v) 


TtT € 8-i/ ^jet 

0 9^^ 9jet 


(4) 


This ultrarelativistic, uniform jet picture is often used in theoret¬ 
ical works about GRBs; 

(iii) relaxing the ultrarelativistic assumption, one must in prin¬ 
ciple take into account the contribution to the observed flux coming 
from the whole emitting volume of the jet (the result of such cal¬ 
culation for the uniform jet is usually dubbed “off-axis jet model”, 
Yamazaki et al. 1 2003|l;|Eichler & Levinson| ( |2004^ ; [Ghisellini| 
2006^; Donaghyl j2006^). For the uniform jet, the result- 




ing apparent structure Eiso(^v) has been computed numerically by 
many authors and it differs from Eq.|^in that the transition from 
the “on-axis” (9^ < 8jet) to the “off-axis” (9^ > 8jet) regime 
is obviously smoother, and a non-zero energy is received from the 
observer even at large viewing angles, since the radiation is not 100 
per cent forward-beamed (red solid line in Fig.|^. 


3.2 A general formula for the apparent structure of a jet 

In the appendix we derive a formula to calculate the apparent struc¬ 
ture Eiao{9v) of ajet with a given (axisymmetric) intrinsic structure 
{e(0),r(0)}. It is valid under the assumptions that the emission 
comes from a geometrically and optically thin volume whose sur¬ 
face does not change significantly during the emission. According 
to our derivation, such apparent structure is given by 



6^9, tf), 9,,) 

m 


e{9) dn 


(5) 


where 9^ is the angle between the line of sight and the jet axis, 
and S is the relativistic Doppler factor. Here Eiao is understood as 
Ttt d|,/(l + 2 ) times the bolometric fluence measured at the Earth 
(dc is the luminosity distance). A formula to calculate the observed 
time integrated spectrum under the same set of assumptions is also 
derived in the appendix (Eq.|33|l. It reads 


E{u, 9^) 


1 + 2 f S^{9,(j),9-,) f{x,a) 

ittdl J r(0) v'Ja 


e{9) dQ 


( 6 ) 


where we have set 2 ; = (1 + z)vl{S vq) for neatness. Here f{x, a) 
is a dimensionless function which defines the comoving spectral 
shape, which can depend on an array a of parameters (see the Ap¬ 
pendix for more details on its definition), izg is some typical fre¬ 
quency of the comoving spectrum, and 


roo 

fs= fix, a) dx (7) 

7o 

Formula]^ can be used to compute the isotropic equivalent energy 
in a specific band, by using 



( 8 ) 


3.3 Comparison with previous studies 

As a consistency check, we test our approach assuming a uniform 
jet structure (eqs. |2| and [3]l and compare the results with the off- 
axis models of [Yamazakiet al.|P003| Y03 hereafter), [ Eichl er ^ 
ILevinsoii] j2004| E04 hereafter) and |Ghisellini et al.| j2006| G06 
hereafter). 
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Figure 3. Apparent structures (Eq.[^ of three uniform jets with F = 100 
and ^jet = 15°, 20° and 30° respectively (solid colored lines). The results 
of |Yamazaki et al.1j2003^ are shown for comparison (solid and dotted thin 
black lines). Our curves are normalized with the same prescriptions as in 
[Yamazaki et al.|f2003t . 


(i) the comparison with Y03 is obtained by using Eq.[^ assum¬ 
ing the same redshift, comoving spectral shape, normalization, and 
Lorentz factor as in Y03. In Fig. we show our results together 
with those of Y03. Apparently, the model used in Y03 (thin black 
lines) slightly underestimates the off-axis isotropic equivalent en¬ 
ergy with respect to ours (colored solid lines); 

(ii) the comparison with E04 is straightforward: the integrand 
of equation 3 of their work, which is used to calculate i?iso(Sv), is 
proportional to <5®, and the same holds for our Eq.|^in the uniform 
jet case. The resulting apparent structures are thus the same up to a 
multiplicative constant; 

(iii) similarly, the integrand of equation 2 of G06, which is used 
to calculate the observed time-integrated^ectrum, is proportional 
to 5^. Again, the same holds for our Eq.lmin the uniform jet case. 
Integration of either equation over all frequencies gives an addi¬ 
tional 5 factor so that, as in the previous case, the resulting apparent 
structures are the same up to a multiplicative constant. 

We conclude that our model is reasonably consistent with previous 
studies on the uniform jet model, and it has the advantage that it 
can be applied to the SI case using the definition of e(0) commonly 
found in the literature. 


3.4 Application to power law and Gaussian jet models 


Here we want to show how the apparent structure of a SJ depends 
on the energy profile e{6) and especially on the Lorentz factor 
profile r(0): the latter has been assumed constant l Graziani et al.| 
2006|> or its role has been deemed secondary jRossi et al.||2002| 


Zhang & Meszaros|2002| ) in many preceding works. 

Fig.|^shows the computed apparent structures of three power 
law jet models (upper panel) and three Gaussian jet models (lower 
panel). The intrinsic structures are defined following |Kumar & Gra-| 
|not| ( |2003t as 


e(0) = 


{e/e. 


e^e, 
)■“ e>e. 


r, e 4 , e, 

1 + (r, -1) {e/e,)-^ e>e. 


( 9 ) 



1 10 100 


0v [deg] 


Figure 4. Apparent structures, according to Eq.|^ of three power law (upper 
panel) and three Gaussian (lower panel) jet models. All power law models 
have Oc = 2°, Fc = 400, a = 2 and different values of the b parameter, 
listed in the legend. The Gaussian models have Oc = 2° and different values 
of Fc, listed. The corresponding intrinsic energy structures (understood as 
47re(0v)) are shown (black dashed lines). The apparent structures in the 
Gaussian case decrease as E[so ~ 9^^ at large viewing angles, regardless 
of the core Lorentz factor Fc and the jet typical angular size Oc- 


e{e) = 

r( 6 ») = 1 -b (Fc - 1) 


( 10 ) 


for the power law and Gaussian jet model respectively. In both 
cases, the 0c parameter represents the typical angular scale of the 
intrinsic structure, i. e. the angle within which most of the jet energy 
is contained. Inspection of Fig. shows that the more the Lorentz 
factor varies, the less the apparent structure mimics the underlying 
intrinsic structure. The Gaussian jet model, in particular, displays 
an apparent structure which is quite well described by a power law 
with a slope around —3, plus a roughly constant core. 


3.5 Reformulation of the Gaussian intrinsic structnre 


In the Gaussian case, the slope of the power law tail of the apparent 
structure at large viewing angles (see Fig. is almost unaffected 
by changes in the two parameters 0c and Fc, and it remains between 
—3 and —4 for reasonable values of these parameters. A different 
slope can be achieved by modifying Eq.[^as follows 


e(0) = 

F(0) = 1 -b (Fc - 1) 


( 11 ) 
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Figure 5. Apparent structures of three Gaussian jets with different values 
of /6c (reported near each line), together with the apparent structure of a 
uniform jet (blue solid line). The con'esponding intrinsic structures are also 
shown (dashed lines). The Gaussian jets have 6c = 3° and Fc = 400; the 
uniform jet has 6jet = 3° and F = 400. 


where Or is a new parameter that allows for the Lorentz factor 
structure and intrinsic energy structure to fall off over different an¬ 
gular scales. In principle 6c might differ from 9r for the following 
reason: t{9) is related to the energy density u = pc? -\- p + ub 
(where ub is the magnetic energy density) of the jet according to 




AtvR'^{9)A{9) 




u{0) 


( 12 ) 


where A{9) is the width of the emitting volume and R{9) defines 
its surface. The energy density u is related to the comoving one by 
u{9) = r^(6)u'(6). Let us take the simplest picture as an example: 
(a) let the emitting volume be a portion of a spherical shell, with 
fixed width A and radius R\ (b) let the efficiency p be the same 
at all angles. One then gets e{9) oc u{9) = T^(9)u'{9). If u' is 
constant, this implies 9r = v/26c- The efficiency, geometry and 
energy density all play a role in determining the ratio 6r/6c. This 
ratio is the main parameter affectiM the slope of the power law tail 
of the apparent structur^ Figurej^shows the apparent structure of 
the Gaussian jet for three values of 9r/9c, along with the uniform 
jet for comparison. 


4 THE GAUSSIAN SJ AS A QUASI-UNIVERSAL JET 
MODEL 

4.1 Luminosity function 

The luminosity function of GRBs clearly depends upon their ap¬ 
parent structure rather than on the intrinsic one. We (jPescalli et al.| 
|2015| > recently showed that a quasi-universal jet with a power law 
apparent structure is consistent with the observed luminosity func¬ 
tion (LF) of Long GRBs: according to our analysis above, the cor¬ 
responding intrinsic structure might then be Gaussian. In order to 
test the possibility that a Gaussian quasi-universal SJ is compati¬ 
ble with the Long GRB LF, we first need a way to relate Eiso and 

^ Let us remark that to get a Gaussian apparent structure, as assumed e. g. 
in |Zhang et al.|j2C)04j and |Dai & Zhang|j2005^ , one needs 6r 3> 6c, i. e. 
the Lorentz factor should remain very high while the energy per unit solid 
angle falls off exponentially. 


Figure 6. The light blue line represents the luminosity function of our fidu¬ 
cial model (see Table [^. There are two shaded areas obtained by vary¬ 
ing one parameter and leaving the others fixed: the gray one refers to 
f 6p/6c ^ \/2, while the pink one refers to cr = ^ + rrj^og ^ 

between 0.35 dex and 0.78 dex. The dashed gray line (visible on the bot¬ 
tom right comer) is the LF for 6c = 5°, while the dashed light blue line 
is the LF for 6r/6c = 0.5, both with all other parameters fixed. The data 
points are the s ame as i n|Pescalli et al.H2 015|, and are partly based on pre¬ 
vious works by jSoderberg et al.H2Q06) and [Wanderman & PiranH2010^ . 
Red points are lower limits. 


the isotropic equivalent luminosity Liao- The simplest approach is 
to define a rest frame duration of the hurst t and to assume a tri¬ 
angular shape for the prompt emission light curve. One then has 


2 X Eiso/t 


(13) 


To define our candidate quasi-universal Gaussian SJ we need 
a set of typical values of the model parameters. Here is an educated 
guess, based on heuristic arguments: 


(i) we define the typical rest frame duration {t) as the mode 
of the observed prompt emission time distribution (Tgo) ~ 70s 


i|Sakamoto et al. 

2011) divided by the average Long GRB redshift 

(l-Fz) « 3.14 

Hjorth et al. 2012 i, obtaining (t) ~ 22s; 


(ii) the Long GRB luminosity function ([Wanderman & Piran| 


2010 1 breaks around L* « 3 x 10° erg/s: in the SJ picture, this 


SEl 


luminosity corresponds to the typical GRB seen on-axis. By Eq.|13| 
such GRB has an on-axis isotropic equivalent energy Ei^oiOv ~ 
0) « 3 X 10®® erg, or equivalently (cc) ~ 3 x 10®®erg/47r « 
2.4 X 10®® erg/sr. Since the highest measured Eiso so far is 
proximately 5 x 10®"^ erg (e. g. GRB 080916C - Ghisellini e 
|2010| >, the Cc parameter requires some dispersion to accommodate 
the observations; 

(iii) the total emitted energ)[^(during the prompt) is 


c7r/2 

Ej = 2 tv e(6) s'mOdO ~ ittcOZ (14) 

Jo 

According to |Kumar & Smoot| ( |2014^ , a typical jet employs around 
10®® erg to break out the envelope of the star in the collapsar sce¬ 
nario: requiring the remaining energy Ljet = E^jp to be at least of 
the same order, assuming an efficiency p = 0.2 we obtain a lower 


Some authors assume two equal, oppositely directed jets. To avoid con¬ 
fusion, we stress that we refer here to a single jet. 
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Table 1. Set of parameters defining our quasi-universal Gaussian SJ. 


Parameter 

Value 

Comment 

Limits“ 

IttCc 

3 X 10®® erg 

Needed to match the break in the observed LF. 

CTiog Cc $ 1 

t 

22s 

Mode of the observed Tgo divided by average redshift {1 -\- z). 

0.35 dex < (Tiog $0.78 dex 

Fc 

800 

Highest r inferred from the onset of an afterglow light curve. 

Fc > 100 

Sc 

3° 

Gives a total energy of the same order of the break out energy. 

3° < 0c < 5° 

er/0c 

1 

Reasonable if the density of the jet core is less than that of the wings. 

0.5 < 0r/0c < V2 


“ all the limits are discussed within the text. 


limit on the jet angular scale 

dc > 3° 


(15) 


A jet with an aperture smaller than this must have lost more than 
half of its initial energy in the excavation of its channel through the 
star envelope; 

(iv) some mixing is likely to occur between the jet borders and 
the stellar envelope ( Rossi et al.|2002| >, and indeed simulations in¬ 
dicate (e.g. |Morsony et al.| i 2007| )) that the jet plasma density in¬ 
creases with t he distance from the axis. In the simplest case out¬ 
lined in (3.5 this suggests a ratio 0r/0c < \/2. We thus take 


dr /0c = 1 for simplicity; 

(v) the exact value of Tc has little effect on the apparent struc¬ 
ture of the Gaussian jet (as long as it is > 100), so it is a secondary 
parameter for what concerns the LF. By the way, let us remark that 
hydrodynamic simulations by |Kumar & Gran^j2003^ suggest that 
for a Gaussian SJ the afterglow onset time is soonest for 0v = 0, 
in which case it is the same as that of a uniform jet with F = Fc. 
Within the SJ picture, this indicates that the highest Lorentz factors 
inferred so far (those of GRBs detected in the GeV energy range by 
Fermi-LAT, see [Ghirlanda et aL]p012[ >) are a measure of the core 
Lorentz factor Fc of the underlying jets. We therefore give this pa¬ 
rameter the rather high value Fc = 800. 

The set of typical parameter values that we induced from these ar¬ 
guments is given in Table [T] As stated in point (ii) above, some 
dispersion in Cc around its typical value is necessary to match the 
observations. Similarly, observations show that the rest frame emis¬ 
sion time t ~ Tgo/l -I- 2 is certainly not universal (e.g. Sakamoto 


jet al.|20lT| l: in the next section we will discuss how to handle these 
two parameters and their dispersion. 

Before moving on, let us note that we are assuming no evolu¬ 
tion of the typical values with redshift: this might well be a rough 
approximation, since the overall progenitor properties may vary 
with redshift. Nevertheless, to keep the discussion as simple as pos¬ 
sible, we neglect this aspect and assume that the quasi-universal jet 
is the same at all times in the past. 


4.1.1 Dispersion on emission time and jet total energy 

As a starting point, we assume a lognormal distribution for both 
the intrinsic duration t and the core energy parameter tc with a 
total dispersion parameter criog ~ 

persion on the luminosity is affected the same way by the disper¬ 
sion on t and on Cc because of Eq.[^. Figure]^ shows the LF of 
the model assuming (Tiogiigc = 0.56 dex (light blue solid line), 
together with the data used in [Pescalli et al.j ( [2015[ l. This is the 
value of (Tiog for which we find the best agreement between 
the theoretical LF and the observed one. Varying aiog between 
0.35 dex and 0.78 dex one still obtains LFs (pink shaded area in 


Fig.|§ that lie within the error bars of the data points, so we take 
these values as estimates of the limits of this paramete^ 

In a realistic model, one might expect the core energy per unit 
solid angle tc to correlate with the emission time t, so a more rig¬ 
orous approach is to interpret the parameter as the resid¬ 

ual variance of log Cc with respect to its (possible) correlation with 
log t. We have two limiting cases: 


(i) in the “worst” case, a linear relation holds between Cc and t, 
so that the dispersion on t gives no contribution to the dispersion on 
Liao. Let us take the standard deviation of log Tgo (for Long GRBs 
observed by Swift), which is around 0.57 dex (jSakamoto et al.| 
[20TT] |, as an estimate of criog t. The dispersion parameter criog 
(and consequently the dispersion on the logarithm of the jet total 
energy) required in this case to reproduce the LF of Fig. is then 
~ 1.13 dex; 

(ii) at the other end, if t and were independent, a 0.56 dex 
dispersion on log! only would be sufficient to reproduce the LF: 
in other words, a single universal value of the jet energy would be 
consistent with the LF (but not with the observed Liao distribution, 
as noted in the preceding section). 


As a result. Fig. suggests that a quasi-universal value (tc) ~ 
2.4 X 10®^ erg/sr (which corresponds to Ljet = 10®® erg if Oc = 
3° and p — 0.2) with a dispersion of less than 1 dex (but not much 
less) is compatible with the considered observational constraints. 
This goes well along with the fact that the progenitor (a former 
Wolf-Rayet star) is expected to possess rather standard properties. 

At this point, a remark is necessary: we assumed that the dis¬ 
tribution of t is independent from the luminosity, which might be 


questionable. Indeed, some authors argue (Daigne & Mochkovitch 


2007l[Virgili et al.[2009[|Bromberg et al.|201 l[|Nakar[2015[|Bamiol 

Duran et al.[2015[ ) that low luminosity, very long GRBs might rep¬ 
resent a distinct population, possibly originating from a different 
progenitor. One of such GRBs is included in the low luminosity 
bin of Fig.|^ namely GRB060218 (the only other burst in the bin 
is GRB980425, which lasted ~ 40s). Based on such a distinction, 
one may suppose that there is an anti-correlation between lumi¬ 
nosity and duration in the overall population. On the other hand, 
though, there is another subclass of GRBs (the so called “ultra- 
long” GRBs, see [Piro et al.|2014{ [Levan et al.|20l4) which have 
durations around 10^ s, but are not underluminous: indeed, no clear 
correlation exists between luminosity and duration within today’s 
samples. We thus conclude that, for our simple model, the assump¬ 
tion of a distribution of emission times which does not depend on 
the burst luminosity is acceptable. 


® Let us remark that the large uncertainties in the obseived LF make it not 
very constraining: as a matter of fact, the reduced Chi-square d of the model 
in Fig.|^is formally yf ~ 10“^. Indeed, we have shown in [Pescalli et al.j 
|2015| that several jet models can reproduce the observed LF to date. 
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4.1.2 Limits on the ratio 9r/O^ 

While the dispersion criog affects the high luminosity end of the 
LF in Fig. the value of the ratio Or/Oc has its influence on the 
low luminosity end. As explained in §3.5[ different values of the 
ratio yield different slopes of the apparent structure. The steeper 
the slope (the higher Or / 0^), the fainter the jet when it is seen un¬ 
der a 90 degrees viewing angle: this implies a lower limit on Or / 0^ 
if we require the predicted LF to extend down to the lowest ob¬ 
served luminosities. For our model, this lower limit is Or jO^ ~ 0.5 
(light blue dashed line in Fig.[^; increasing OrjO^ from 1 to \f2 
one obtains the family of LFs spanning the gray shaded area in 
Fig. The highest value of Or/O^ consistent with the lower limits 
on the rate of intermediate luminosity GRBs (red points in Fig.|^ 
is Or/Oc ~ x/2. 


3.5 

3 



1 



Log(Eiso) [»g] 


4.1.3 Limits on the angular scale Oc 

The width of the angular scale Oc impacts mainly on the high lu¬ 
minosity end of the LF. The wider Oc, the higher the probability of 
observing the jet within the core (i. e. ^ Oc), implying a higher 
rate of bursts with high observed luminosity. The dark grey dashed 
line in Fig. shows how the LF would change if Oc = 5° was as¬ 
sumed. A wider Oc would lead to higher rates than those observed 
at high luminosities. 

As explained in ^4.1[ the lower limit Oc ~ 3° is given by 
the requirement that the jet total energy Ej^t = E^/ri of the quasi- 
universal jet is at least 10®'^ erg, i. e. of the same order of the energy 
previously spent by the jet to excavate its way through the stellar 
envelope ( [Kumar & Smoot|2014| l. This is not a strict requirement, 
since one can also have that most of the jet energy is spent prior 
to the break out (we stress that Ejet is the total jet energy after the 
break out). 


4.2 Consistency with the Amati correlation 

Gamma Ray Bursts show a strong correlation j Amati et al.|2002^ 
[Lamb et al.|2005[|Lloyd et al.|2000^ between the peak of their oEi, 
spectrum (i?peak) and the isotropic equivalent energy Eiso, roughly 

Fipeak OC E^/j (16) 

A similar correlation involves the isotropic equivalent luminosity 


Liao ( jYonetoku et al.|[2004[ (. These correlations have been exten¬ 
sively studied in the past years for (a) their possible physical im- 


plications on the GRB emission process (e.g. Rees & Meszaros 

(2005 

; Toma et al. 

(|2005|i; |Barbiellini et al. (|2006|i; |Ryde et al. 

(2006 

; Thompson[ 

2006f; [Giannios & Spruit (2007i; Thompson 

et al.[ 

2007J) and on the GRB jet structure ([Yamazaki et al.|2004| 

Eichler & Levinson|2004| |Lamb et al.|2005l|Levinson & Eichler 

2005|, (b) the possibility to use them to standardize GRB ener- 


getics and constrain the cosmological parameters ( [Ghirlanda et al.j 
|2004[|Amati et al.|2008|l. Besides, these correlations stimulated an 


intense debate on the possible impact of selection effects ([Nakar & 


Piran|2005|[Band & Preece|20051 [Butler et al.|2007|[2009||Shah- 

moradi & Nemiroffl2009 

|Kocevski||2012 

Ghirlanda et al.|2005| 

Bosnjak et al.|2008[ Ghir 

anda et al.|2008 

Nava et al.|2008[[Am- 

ati et al. [20091 [Krimm et al.|20091|Ghirlanc 

a et al.|2012||Heussaft 


et al.|2013|l. Despite the wealth of papers, the spectral-energy cor¬ 


relations are still a hot subject in the field and no consensus on their 
physical interpretation has been reached yet. 

It has been proposed recently that the Amati correlation might 


Figure 7. The thick solid colored line represents a sequence of Liao(®v) 
and Speak (Ov) for a Gaussian SJ with the parameters given in Tab^ The 
color coding of the curve accounts for the viewing angles according to the 
scale defined on the right. Data points from the Swift BAT complete GRB 
sample [Nava et al.|2012t are shown for comparison. The pink dashed lines 
represent the Amati coiTelation and the 3a dispersion of the data points as 
computed by jNava et al.jjlOT^ . The shaded area represents the portion of 
the plane spanned by the curves {Siso(®v)) Speak(l5v)} with log Cc vary¬ 
ing within itcr (darkest shade), ± 20 - and itScr (lightest shade) respectively 
(we assumed a = 0.5 dex). 


be due to a sequence of bulk Lorentz factors, with more luminous 
GRBs having larger F values, as suggested by the possible clus¬ 
tering of the GRB energetics when transformed in the comoving 
frame ( [Ghirlanda et al.|2073) . 

Here we wish to show that, within the quasi-universal Gaus¬ 
sian SJ defined in the preceding section, the Amati correlation can 
be explained as a viewing angle effect. First, we need to make some 
assumptions on the comoving spectrum emitted by the jet: 


(i) following e. g. [Yamazaki et al.[p003[ l; [Ghisellini et al.[[2006[ l 
we assume a smoothly broken power law shape for the comoving 
spectrum. In the notation of Eq.[^ this amounts to write 


f{x,a,b) 


1 

a;“ -I- I*’ 


(17) 


We choose a = —0.1 and b = 1.3, which are the typical low 
and high energy slopes of observed GRB energy spectra (e. g. [Nava[ 
[et al.[ [201 l[ l) (the results of this work do not change if we take 
different values of a and b within their typical dispersion); 

(ii) we choose the peak frequency of the comoving spectrum to 
he Oq — 1 keV at all points in the jet. This is motivated by the 
finding ( [Ghirlanda et al.[2012^ that the peak of GRB spectra seems 
to cluster in the comoving fram^ To choose a different value for 
u'q amounts to move vertically the colored curve in Fig.|^ 

With the above assumptions, we used Eq. to calculate the ob¬ 
served spectrum and Eq. to compute the corresponding Eiso in 
the [1 keV — 10"* keV] band. The solid colored line in Fig. [^rep¬ 
resents the resulting sequence of Eiao{0v) and i?peak(0v) for our 
Gaussian SJ with the parameters given in Tab. Data points (for 


® One could also take Uq = 0^(0), i.e. give a structure to the peak of 
the comoving spectrum, and this would certainly lead to different results. 
Since we have no argument to prefer such a choice, we limit ourselves to 
the simpler case discussed in the text. 
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comparison) are bursts from the Swift BAT complete sample jSal-| 
|vaterra et al.|2Q12) which is flux limited and 97 per cent complete 
in redshift; the pink dashed lines represent the Amati correlation 
and the Sa dispersion of the data points as computed by |Nava et al.| 
( |2012^ using the same sample. Changing the value of (a dis¬ 
persion of Cc is necessary to reproduce the LF, as discussed in the 
preceding section) amounts here to moving horizontally the colored 
curve: assuming criogt^ = 0.5 dex, we colored the portions of the 
plane spanned by the curve when log Cc varies within ±(j (darkest 
shade), ±2 (t and ±3 (j (lightest shade) respectively. 

We conclude that this model is compatible with the interpre¬ 
tation of the Amati correlation as a sequence of viewing angles, 
with the dispersion of the observed points being due to the intrin¬ 
sic dispersion in the total jet energy and possibly in the peak of the 
comoving spectrum. 


5 DISCUSSION AND CONCLUSIONS 


To some extent, structure is an inevitable feature of GRB jets, since 
a uniform jet with discontinuous edges is clearly unphysical: in¬ 
deed, a Gaussian jet structure has been often described as a more re¬ 
alistic version of the uniform jet structure (e. g. [Zhang & Meszaros| 
( |2002 ^) and it is expected to reproduce most of the features of the 
latter when the viewing angle is small enough ([Kumar & Granotj 
|2003| (. 

In this view, our Gaussian SJ model is to be understood as a 
proxy for a jet in which both the emissivity and the Lorentz factor 
decrease rapidly away from the axis. Our results then indicate that 
the simplified “ultrarelativistic uniform (top-hat)” jet model, while 
very useful, likely predicts too little off-axis energy emission with 
respect to any more realistic counterpart. This might be a clue to 
the still missing observation of “orphan afterglows”: the jet can be 
narrow (which is a necessary condition, in the collapsar scenario, 
for the jet to break out the progenitor star envelope without wasting 
too much energy), and still be visible at viewing angles larger than 
its typical angular dimensioij^ contrary to the uniform jet model 
( [Perna & Loeb|19^ . This affects also the interpretation of the lu¬ 
minosity function of GRBs: the high rate of underluminous events 
and the wide range of observed luminosities (from 10^^ erg/s up to 
10®^ erg/s) are readily explained without invoking the existence of 
different burst populations ( [Pescalli et al.[2015| l. 

How standard can be the jet of GRBs in this picture? The ob¬ 
served break in the luminosity function of Long GRBs at L* ~ 
3 X 10®^ erg/s i Wanderman & Piran|2010 sets a natural scale for 
the luminosity. Within the quasi-universal structured jet hypothesis, 
this luminosity corresponds to the typical jet seen on-axis. The on- 
axis Eiso of the typical jet is thus i?* ~ 3 x 10®® erg (the jet total 
energy Fljet = E-y/r) ~ 2'K6^ec depends on the angular dimen¬ 
sion 8c of the jet, and it is Ejet ~ 10®® erg for dc = 3°). The fact 
that we do observe more energetic GRBs (up to Eiso ~ 5 x 10®'® 
erg) means that some dispersion in the jet energy is necessary to 
account for it. On the other hand, a dispersion of around 0.6 dex in 
the maximum luminosity L* is needed to reproduce the luminosity 
function. As discussed in j [4.1.1[ this limits the dispersion on the 
jet total energy below 1 dex. Figure|^shows that 0.5 dex is enough 
to account for the dispersion in the Amati correlation. We can con¬ 
clude that the jet of GRBs in this picture can be rather standard. 


^ The reduction of the expected number of orphan afterglows is a common 
feature of SJ models, see|Rossi et al.|[2008|. 


with a total energy that lies within a factor of 10 from the typical 
Ejet ~ 10®® erg in most cases. This is exactly what one would 
expect by the association of GRBs with supemovae Ib/c, since the 
latter should have rather standard progenitors. 


The explanation of the Amati correlation as a viewing angle 
effect has been proposed several times in the past, within a vari- 


ety of jet models (e. g. Yamazaki et al. \ 2004|l; 

Eichler & Levinson 

(|2004J; [Lamb et al.j (|2005|l; [Levinson & Eich 

er| ([2005^; [Graziani 

et al.|(200^; recently also in a photospheric 

emission model by 


Lazzati et al.|(|2013ll, where the jet structure is computed through 


relativistic hydrodynamical simulations of the jet break out). The 
result presented in j ]4.2[ shows that this interpretation is possible 
also in our simple model, at least within the discussed assumptions. 
In this view, the dispersion of the observed correlation is due to the 
intrinsic dispersion on the jet total energy (and possibly on the co¬ 
moving peak energy). 

We conclude by stressing (as noted in §4.1[ that in this simpli¬ 
fied model we neglected the possible evolution of the universal jet 
parameters with redshift, which might play an important role in de¬ 
termining the luminosity function and other observational features 
of the GRB population. We also assumed no correlation between 
emission time and luminosity, for the reasons explained in the last 
paragraph of §4.1.1[ 
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6 APPENDIX 

6.1 Derivation of the formula for Eiso {0 ^) 

As a first step, consider an uniformly expanding hemispherical shell which 
emits radiation for a small time interval At, during which the radius does 
not vary appreciably (see Fig. |8]for a sketch of the geometry). For defi¬ 
niteness, we set up two coordinate systems: the origin of the first (call it 
K) is the center of the hemisphere, and its z axis lies on the line connect¬ 
ing this point to the observer. The spherical coordinates of this system will 
be referred to as (r, 0, 0). The second system (A'l) is centered on the ob¬ 
server, and its zi axis coincides with z, but it is oppositely oriented. The 
con'esponding spherical coordinates will be called (ri, 0i). If there is 
no significant absorption, the flux received by an observer at a distance d is 
(neglecting cosmological corrections) 

= — [ d(pi s'mOidOiridri — ri/c) (18) 

Js(t) 

where S (t) is the equal arrival time surface for photons received at time t, 
andju is the specific emissivity. To obtain the corresponding fluence 7^(1^), 
we must integrate over time. This allows us to “rearrange” the emission 
times so that we do not need to bother about the equal arrival time surfaces, 
i. e. 

E{i') = ^ [ dt [ d<pi smOidOirfdri (19) 

JAt Jv{t) 

where V (t) is the emitting volume at time t. Now we introduce some sim¬ 
plifications: 
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(i) the hemisphere is far away from the observer {d ^ R), so that we 
can safely set ri 5^ d, sin 0i ^ 6i and 9i ^ sin 9 K/d\ 

(ii) since we integrate (j)i from 0 to 27r, we can equivalently integrate 
over (p\ 

(hi) we assume that the emitting volume is geometrically thin, i. e. it is 
a hemispherical shell of width Ar. 


Since dri cos 9 = dr, we have 

d2 r ^7r/2 rlir r 

T{y) = —;r / dt j sin9d9 / / ju{F,t)dr (20) 

d^ JAt Jo Jo Jav 


Now we use the relations ji, = j'^, jRybicki & Lightman|l979 


where <5 


is the relativistic Doppler facto^and primed quantities refer to the comoV' 
ing frame, and dr = Sdr' to integrate over Ar' and obtain 


where we assumed (without loss of generality) that the line of sight lies on 
the z — X plan^ 

Finally 

Eiso(0v) = J de) dfi (31) 

Summaiizing, the formula above gives the apparent stnacture of a jet, 
given its intrinsic structure (i. e. e{9) and r(0)), seen under a given viewing 
angle 9v, under the assumption that the emission comes from a thin, trans¬ 
parent volume, whose surface R{9, (p) does not change significantly during 
the emission. 


6.2 Derivation of the formula for T{u,0v) 


d 2 n n'n II nl'K 

Tiy) = —— / dt sin9d9 / (0,0, f) (21) 

JAt Jo Jo 

where , is the comoving specific intensity. Integrating over dt = dt' jb 
and dv = bdv' to get the bolometric fluence, we have 

d2 r27T /’7r/2 

T=— / < I'> At'sin ed9dif> (22) 

d^ Jo Jo 

where < 7' > is the average comoving intensity during the emission time 
If the sphere emits uniformly, < I' > does not depend on 9 and 0: in 
this case the integral is analytic, yielding 


d^ 

By definition we have Eiso = 47r d^E, so we finally get 

= 47r2(l + 13 f {2 - j3)R^ T < l' > At' 


(23) 

(24) 


As stated in section[3^ in the ultrarelativistic limit the isotropic equivalent 
energy is indistinguishable from that of a spherical explosion, which yields 
E[so = 47re by definition, thus in this case we have 

IQn^R^ r < I' > At' = Ane (25) 


so that we can make the identification 

<I' > At' = -^ 

47r i?2r 


(26) 


Since the shell is geometrically thin, the intensity coming from a point 
(0, (p) is due only to the local emitting volume, so in this approximation 
we can think of Eq.[^as a relation between local quantities, namely 

,27, 

We thus substitute this equivalence back into Eq. and multiply it by 
47r d ^, to get 


E- — 

-^ISO - 


^27r p- 

Jo Jo 


b^{9,(P) 

no.<p) 


(p) sin 9 d9 dtp 


(28) 


In words, this equation tells us how to weigh the contribution from each el¬ 
ement of the hemisphere in order to take into account the relativistic effects 
and the local energy density e(9, tp). 

As long as the expansion is purely radial, this equation holds for el¬ 
ements of any surface - in other words, setting R = R{9,(p) does not 
affect the validity of the derivation. In the case of a structured jet we have 
e = €(9), r = r(0) and, if the observer is off-axis, the Doppler factor (5 
must take into account the angle between the line of sight and the velocity 
of each point of the surface. A little geometry allows one to write 


s{e, (j), 0 v) 


with 


1 

r(0Hr^(6») cos a{9, (P, 6>v)] 


(29) 


cos q:(0, 0, ^v) = cos 9 cos 9yf + sin 9 sin 0 sin 0v 


(30) 


® We are implicitly assuming that the emissivity is isotropic in the comov¬ 
ing frame. 


We can also derive the corresponding formula for the time integrated spec¬ 
trum E{v) as a function of the viewing angle 9^. First, let us write the 
comoving specific intensity I', as (7 q/z/q) jiy' j^'ot where 7 q is a nor¬ 
malization constant, v'^ is some preferred frequency and j(y'I v'^^ol) is a 
dimensionless function that defines the shape of the comoving spectrum, 
which depends on an array q of parameters. Let us also call the integral 
of f{x^ a) over all positive x's. Then we rewrite Eq. |26| as 

J' = . -.... 1 -- (32) 

^ Aff^A-TTR^F 

Starting again from Eq. |21| taking into account the above definitions we end 
up with 


E{u,9^) = 


1 + z f 5^(19, 0,6»v) /(a:, a) 


A-wdl 


f 


e{9)dQ. 


r(0) v'^fa 

where the comoving spectral shape f(x,ci) is to be evaluated at x 
z)ul(5v'^). 


(33) 
(1 + 


^ Some authors prefer to set the coordinates so that the line of sight lies on 
the y — z plane: in this case, one would have 

cos ol{9^ 0, 0v) = cos 9 cos 0v + sin 9 cos 0 sin 0v 
















